Arbitrary Dimensional Majorana Dualities and Architectures for Topological Matter 
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Motivated by the prospect of attaining Majorana modes at the ends of nanowires, we analyze 
interacting Majorana systems on general networks and lattices in an arbitrary number of dimensions, 
and derive various universal spin duals. Such general complex Majorana architectures (other than 
those of simple square or other crystalline arrangements) might be of empirical relevance. As these 
systems display low-dimensional symmetries, they are candidates for realizing topological quantum 
order. We prove that (a) these Majorana systems, (b) quantum Ising gauge theories, and (c) 
transverse-field Ising models with annealed bimodal disorder are all dual to one another on general 
graphs. This leads to an interesting connection between heavily disordered annealed Ising systems 
and uniform Ising theories with nearest-neighbor interactions. As any Dirac fermion (including 
electronic) operator can be expressed as a linear combination of two Majorana fermion operators, 
our results further lead to dualities between interacting Dirac fermionic systems on rather general 
lattices and graphs and corresponding spin systems. The spin duals allow us to predict the feasibility 
of various standard transitions as well as spin-glass type behavior in interacting Majorana fermion 
or electronic systems. Several new systems that can be simulated by arrays of Majorana wires 
are further introduced and investigated: (1) the XXZ honeycomb compass model (intermediate 
between the classical Ising model on the honeycomb lattice and Kitaev's honeycomb model), (2) a 
checkerboard lattice realization of the model of Xu and Moore for superconducting (p -I- ip) arrays, 
and a (3) compass type two-flavor Hubbard model with both pairing and hopping terms. By the use 
of dualities, we show that all of these systems lie in the 3D Ising universality class. We discuss how 
the existence of topological orders and bounds on autocorrelation times can be inferred by the use 
of symmetries and also propose to engineer quantum simulators out of these Majorana networks. 

PACS numbers: 05.30.-d, 03.67.Pp, 05.30.Pr, 11.15.-q 



I. INTRODUCTION 

Majorana (contrary to Dirac) fermions are particles 
that constitute their own anti-particles.^ Early quests 
for Majorana fermions centered on neutrinos and fun- 
damental issues in particle physics that have yet to be 
fully settled. If neutrinos were Majorana fermions then 
neutrinoless double f3 decay would be possible and thus 
experimentally observed. More recently, there has been 
a flurry of activity in the study of Majorana fermions 
in candidate condensed matter realizations, ^^^® including 
lattice^^"^^ and other^'^ systems inspired by the prospect 
of topological quantum computing. ^'''^^ In the condensed 
matter arena, Majorana fermions are, of course, not fun- 
damental particles but rather emerge as collective exci- 
tations of the basic electronic constituents. The systems 
discussed in this work form a generalization of a model^^ 
that largely builds and expands on ideas considered by 
Kitaev*'^^'^° including, notably, the feasibility of creat- 
ing Majorana fermions at the endpoints of nanowires. ^^ 
A quadratic fermionic Hamiltonian for electronic hopping 
along a wire in the presence of superconducting pairing 
terms (induced by a proximity effect to bulk supercon- 
ducting grains on which the wire is placed) can be ex- 
pressed as a Majorana Fermi bilinear that may admit free 
unpaired Majorana Fermi modes at the wire endpoints. ^^ 
Kitaev's proposal entailed p-wave superconductors.^ 

More recent and detailed studies suggest simpler 
and more concrete ways in which zero-energy Majo- 



rana modes might explicitly appear at the endpoints of 
nanowires placed close to (conventional s-wave) super- 
conductors. Some of the best known proposals^'^'^*^ en- 
tail semiconductor nanowires (e.g., InAs or InSb^'^) with 
strong depolarizing Rashba spin-orbit coupling that are 
immersed in a magnetic field that leads to a competing 
Zeeman effect. These wires are to be placed close to su- 
perconductors in order to trigger superconducting pair- 
ing terms in the wire. By employing the Bogoliubov-de 
Gennes equation to study the band structure, it was read- 
ily seen how Majorana modes appear when the band gap 
vanishes. ^'^'^° Along another route, it was predicted that 
zero-energy Majorana fermions might appear at an in- 
terface between a superconductor and a ferromagnet.^'^® 
Majorana modes may also appear in time-reversal invari- 
ant s-wave topological superconductors.^^ 

If zero-energy Majorana fermions may indeed be har- 
vested in these or other ways^ then it will be natural 
to consider what transpires in general networks made of 
such nanowires. The possible rich architecture of struc- 
tures constructed out of Majorana wires and/or particu- 
lar junctions may allow for interesting collective phenom- 
ena as well as long sought topological quantum comput- 
ing applications. ^°'^^ Interestingly, as is well appreciated, 
the braiding of (degenerate) Majorana fermions realizes a 
non-Abelian unitary transformation that may prove use- 
ful in quantum computing providing further impetus to 
this problem. In the current work, we consider general 
questions related to Majorana Fermi systems that may 



be constructed from nanowire architectures. 

A central question regarding systems of Majorana 
fermions is concerned with viable topological quantum or- 
ders (TQOs). Disparate (yet inter-related) definitions of 
TQO appear in the literature. One of the most strik- 
ing (and experimentally important) aspects of TQO is 
its immunity to local perturbations or, equivalently, its 
inaccessibility to local probes at both zero and finite 
temperatures.^^ Some of the best studied TQO systems 
are Quantum Hall fluids. ^^ Several lattice models are also 
well known to exhibit TQO, including the spin S = 1/2 
models introduced by Kitacv.^^'^° In the context of the 
Majorana lattice systems (and general networks) that we 
investigate here, one currently used approach for assess- 
ing the presence of TQO^^ is observing whether a for- 
tuitous match occurs, in perturbation theory, between 
(a) the studied nanowire systems with (b) Hamiltonians 
of lattice systems known to exhibit TQO. While such 
an analysis is highly insightful, it may be hampered by 
the limited number of lattice systems (and more general 
networks) that have already been established to exhibit 
TQO. 

In this work we suggest a different method for 
constructing Majorana system architectures displaying 
TQO. This approach does not require us to work to- 
wards an already examined lattice system that is known 
to exhibit TQO. Instead, our recipe invokes direct conse- 
quences of quantum invariances. Symmetries can man- 
date and protect the appearance of TQO^"*^ via a gen- 
eralization of Elitzur's theorem. ^^'^^ Specifically, when- 
ever d- dimensional gauge-like symmetries'^ are present 
(most importantly, discrete d = 1 or continuous d = 1, 2 
symmetries), finite temperature TQO may be mandated. 
Zero-temperature TQO states protected by symmetry- 
based selection rules can be further constructed. A sym- 
metry is termed a d-dimensional gauge-like symmetry if 
it involves operators/fields that reside in a d-dimensional 
volume. ^^^^^ The use of symmetries offers a direct route 
for establishing TQO that docs not rely on particular 
known models as a crutch for establishing its presence. 

To illustrate the basic premise as it may be applied to 
architectures with Majorana fermions, we will advance 
and study a generalization of a model introduced in Rcf. 
19 to describe a square lattice array of Josephson-coupled 
nanowires on superconducting grains. A schematic of the 
array studied in Ref. 19 is presented in Fig. 1. As we 
will elaborate on in Section II, our general-dimensional 
extension of this Hamiltonian is given by Eqs. (5), (6), 
and (8) with cu, (i = 1,2) denoting Majorana operators 
(satisfying the standard Majorana algebra of Eq. (4)) 
associated with nanowire endpoints. Within the gener- 
alized scheme, these nanowires are placed on supercon- 
ducting islands that occupy the vertices r of a general 
(even-coordinated) network, with links I connecting the 
islands. The ends of the nanowires are placed so that 
each link I connects two Majorana fermions cn,cj2 from 
different wires. Each link carries an arbitrary but fixed 
orientation, just for the purpose of labelling the Majo- 



ranas on it: As one traverses a link in the specified di- 
rection, cii comes before Cj2 (see Fig. 1). 




FIG. 1. A decorated square lattice (with unit vectors ei and 
62) in wliicli each site is replaced by a tilted square (repre- 
senting a superconducting grain at site r). Two nanowires 
(solid blue diagonal lines) are placed on each grain. The 
grains are coupled to each other via Josephson couplings. 
A local (gauge) symmetry operator of the model is Gp = 
{iciiici^2){ici^ici^2)iicigicig2){ici2ici^2), where P defines the 
minimal closed loop. See text. 



For example, in Fig. 1, two parallel nanowires are 
placed on each superconducting grain. These grains are 
placed on the sites r of a square lattice matrix. The two 
nanowires on each grain yield four Majorana fermionic 
degrees of freedom, placed on the edges of the oriented 
links of the lattice. The Majorana fermions on different 
superconducting grains, sharing a link, are coupled to 
each other by Josephson junctions. Prior to introducing 
the Josephson couplings, each grain is shunted to main- 
tain a fixed superconducting phase and is capacitively 
coupled to a ground plate. Consequently, there are large 
fluctuations in the electron number operator. However, 
the electron number parity is conserved. The sum of the 
two dominant effects: (i) intcr-grain Josephson couplings 
and (ii) intra-grain constraints on the electron-number 
parity, complemented by exponentially small capacitive 
energies, leads to a simple effective Hamiltonian. The 
intra-grain constraint on electron number (even/odd) 
parity is more dominant than inter-grain effects. The 
parity operator is Vr = (—1)"'' with n^. the total num- 
ber of electrons on grain r. This electron number parity 
can be of paramount importance in interacting Majo- 
rana systems. ^^'^^ In grains having two nanowires each, 
the electronic parity operator is quartic in the Majorana 
fermions; it is just the ordered product of the four Majo- 
rana fermions at the endpoints of the nanowires on top 



of the grain at site r, 

"Pr = Q11Q21Q32Q42, r e li,l2,h,h (1) 

(we write r E I to indicate that r is one of the two 
endpoints oi I). This gives rise to a term in the effective 
Haniiltonian of the form^^ 



Ho = -hJ2'Pr, 



(2) 



with the sum taken over ah grains, whose total number 
is Nr- This term is augmented by Josephson couphngs 
across inter-grain links I, leading to a Majorana Fermi 
bilinear term involving the coupled pair of Majoranas 

{{C11,C12)}, 



h,=-jj: 



lCllCl2- 



(3) 



Fcrmionic parity effects are more dominant than Joseph- 
son coupling {h 3> J) effects. Invoking perturbation 
theory, for small (J/h), it was found^^ that, to lowest 
non-trivial order, the resultant effective Hamiltonian was 
identical to that of Kitaev's toric code model, ^° thus es- 
tablishing that such a system may support TQO. Unfor- 
tunately, for (J/h) ^ 1; spectral gap is small and the 
system is more susceptible to thermal fluctuations and 
noise. A Jordan- Wigner transformation was invoked^^ 
to illustrate that these results survive for finite {J/h). 

In this article, we will outline a general procedure for 
the design of different architectures of nanowires on su- 
perconducting grains that support TQO. As alluded to 
above, our considerations will not be limited to the use 
of perturbation theory but will rather rely on the use 
of symmetries and exact generalized dualities associated 
with these granular and other systems defined on general 
networks. We will further invoke a general framework for 
dualities that does not require the incorporation of known 
explicit representations of a spin in terms of Majorana 
fermions nor Jordan- Wigner transformations that have 
been invoked in earlier works. ^^'^^'^^ The bond-algebraic 
approach, ^^'^^"'^'^ that we employ to study general exact 
dualities and fermionization,'^-'^''^^ allows for the deriva- 
tion of earlier known dualities as well as a plethora of 
many new others for rather general networks (or graphs) 
in arbitrary dimensions and boundary conditions. It is 
important to note, as wc will return to explicitly later, 
that as Dirac fermions can be expressed as a linear com- 
bination of two Majorana fermions, our mappings lead 
to dualities between standard (non-Majorana) fcrmionic 
systems and spin systems on arbitrary graphs in gen- 
eral dimensions. These afford non-trivial examples of 
fermionization in more that one dimension. 

Among several exact dualities that we report here we 
note, in particular, the following: 

• A duality, in any dimension, between the Majorana 
fermion system corresponding to an arbitrary net- 
work of nanowires on superconducting grains and 
quantum Ising gauge theories. 



• A gauge-reducing and emergent dualities in ar- 
bitrary number of dimensions between granular 
Majorana Fermi systems on an arbitrary network 
and transverse-field Ising models with annealed ex- 
change couplings. In two dimensions, this duality, 
along with the first one listed above, indicates that 
an annealed average over a random exchange may 
leave the system identical to a uniform transverse- 
field Ising model. 

• A further duality between a particular Majorana 
fermion architecture and a nearest-neighbor quan- 
tum spin S — 1/2 model which, in some sense, is 
intermediate between an Ising model on a honey- 
comb lattice and the Kitaev honeycomb model. ^^ 
We term this system the "XXZ honeycomb com- 
pass model" . 

As one of the key issues that wc wish to address 
concerns viable TQO, boundary conditions may be of 
paramount importance. Boundary conditions are inher- 
ently related to the character (and, on highly connected 
systems, to the number) of independent d-dimensional 
gauge-like symmetries. Imposing periodic or other 
boundary conditions on a system can lead to vexing prob- 
lems in traditional approaches to dualities and fermion- 
ization. By using bond-algebras, we can circumvent these 
obstacles and construct exact dualities for both infinite 
systems and for finite systems endowed with arbitrary 
boundary conditions. Other formidable barricades, such 
as the use of non-local string transformations, can be 
overcome as well within the bond-algebraic approach to 
dualities. "^^ The validity of any duality mapping can, of 
course, be checked numerically by establishing that the 
spectra of the two purported dual finite systems indeed 
coincide. The matching of the spectra serves as a defini- 
tive test since dualities are (up to global redundancies) 
unitary transformations'^^ that preserve the spectrum of 
the system. 



II. NETWORKS OF SUPERCONDUCTING 
GRAINS AND NANOWIRES 

In the Introduction, we succinctly reviewed the effec- 
tive Hamiltonian for the square lattice array, ^^ depicted 
in Fig. 1, of Josephson-coupled granular superconductors 
carrying each two nanowires. This architecture serves as 
a useful case of study. There is more to life, however, 
than square lattice arrays (although we will return to 
these later on in this work). We consider next rather 
general architectures in which each node r (supercon- 
ducting grain) has an even number of nearest neighbors 
to which it is linked by Josepshon coupling, see Fig. 2. 
These general networks include, of course, any two di- 
mensional lattice of even coordination, e.g, those of Figs. 
I, and 3, as special cases. 

The architectures that we consider are realized by plac- 
ing at each vertex r of a graph-theoretical network a 




FIG. 2. A general network of superconducting grains with an 
even coordination number of each vertex. The local coordina- 
tion number q^ of any superconducting grain centered about 
site r is equal to the number of endpoints of all nanowires 
that are placed on that superconducting grain. The domi- 
nant Josephson tunneling paths between inter-grain nanowire 
endpoints are highlighted by solid lines. Shown here is a two- 
dimensional projection of the network. 



finite-size superconducting grain. On each of these grains 
there are Zj. nanowires. These nanowires provide 2zr Ma- 
jorana fermions, one for each wire's endpoint. Inter-grain 
Josephson tunneUng is represented by a link involving 
Majoranas coming from different wires on different is- 
lands. We place the nanowires on every grain in the 
network so that each endpoint of a nanowire is near the 
endpoint of another nanowire on a neighboring grain, to 
maximize Josephson tunneling. Thus, the coordination 
number q^. of grain r in these graphs is g^ = 22^.'^^ The 
general situation is depicted in Fig. 2. 




• an intra-grain charging energy hr associated to 
each island at site r. 

In a general, spatially non-uniform, network the spatial 
distribution of couplings Ji and charging energies hr need 
not be constant. 

The algebra of Majorana fermions is defined by the 
following relations: 



{cii,ci'i'} = ^Si^i'St^i 



Cli 



(4) 



With all of the above preliminaries in tow,'^^ we are now 
ready to present the effective Haniiltonian for the systems 
under consideration, 



Hm = -i/^ Ji'^'ii^i'i ^ zJ ^rVr 



where 



I J* — % ^llil ^^2^2 



r eh,-- ,L,, 



(5) 



(6) 



is the product of all Majorana fermion operators asso- 
ciated with the superconducting grain at site r, ordered 
in sonic definite but arbitrary fashion (differing ordcrings 
produce the same operator up to a sign)."^^ 

The index im can be either i^ = 1 or i,„ = 2, depend- 
ing on the particular orientation that has been assigned 
to the links in the network. More precisely, im = 1 if Im 
points away from r, and i^ = 2 if Im points into r. The 
factor i^^^ is introduced to render Vr self-adjoint. Since 



(cjin 



■Ci„ 



Jt ^ (_l),.(,.-l)/2,^^^^ 



■Ci„ 



(7) 



and Qr = 2zr, we set the integer Zr2 to be the number of 
nanowires counted modulo 2, 



if Zr is even 

1 if Zr is odd 



(8) 



As we remarked earlier, the operators Vr are related to 
the operators rir counting the total number of electrons 
on the grain r as 



Vr = (-l)"^ 



(9) 



thus measuring the parity of the number of electrons at 
site r. Haniiltonian (5) constitutes an arbitrary dimen- 
sional generalization of the sum of the two terms in Eqs. 
(2, 3). In the following, we call the operators {icnQ2} 
and {Vr} the bonds of the Hamiltonian Hm-^^'^^ 



FIG. 3. A triangular network of superconducting grains 
(hexagons) on each of which we place three nanowires. 

The basic inter and intra-island interactions have dif- 
ferent origins. For ease of reference, we reiterate these 
below for arbitrary networks: 

• there is a Josephson coupling Ji associated with 
each inter-grain link I of the network connecting 
different superconducting grains, and 



III. SYMMETRIES AND TOPOLOGICAL 
QUANTUM ORDER 

For the particular case of the square lattice {D = 2), 
the interacting Majorana Hamiltonian Hm with periodic 
(toroidal) boundary conditions was found to exhibit 0- 
dimcnsional local, d ~ 1-dimcnsional gauge-like, and 2- 
dimcnsional global symmetries. ^^ These symmetries, in- 
herently tied to TQO^"* and dimensional reduction, ^^"^^ 



also appear in the more general network renditions of the 
granular system just described in the previous section. 
They are also manifest for the interacting Majorana sys- 
tems embedded in any spatial dimension D > 2 when 
different boundary conditions are imposed. '^^ 

Global Symmetry: 
The Hamiltonian Hm of Eq. (5) displays a global symme- 
try Q, given by the product of all the Majorana Fermion 
operators in the system. We can write Q in terms of 
bonds as 



Q = l[Vr, 



(10) 



since each Majorana is contributed by some island. The 
order of the bonds in Q is not an issue, since 



[Pr,rr']^0, 



(11) 



for any pair of sites r,r'. The conserved charge Q repre- 
sents a Z2 symmetry of the system, 



Q2 = l. 



(12) 



Beyond this global symmetry, the system of Eq. (5) 
exhibits independent symmetries that operate on finer, 
lower-dimensional regions of the network. Of particular 
importance to TQO are d = 1 and d ~ O-dimensional 
symmetries, and so we turn to these next. 

d = 1 symmetries: 
The d = 1 dimensional symmetry operators of the Majo- 
rana system are given by 



Qe = Y[{iciici2), {Qif = 1, 



(13) 



lai 



where ^ is a continuous contour, finite or infinite and 
open or closed depending on boundary conditions, en- 
tirely composed of links. That these non-local operators 
are symmetries is readily seen once it is noted that (a) 
each of the terms (or bonds) in the summand of Eq. (5) 
defining Hm involves products of an even number of Ma- 
jorana fermions and (b) by the second of Eqs. (4), effect- 
ing an even number of permutations of Majorana fermion 
operators in a product incurs no sign change. For exam- 
ple, for a network of linear dimension L along a Cartesian 
axis, the contour i spans 0{L^) sites and is thus a d = 1 
dimensional object. This is the origin of the name d = 1 
symmetries. Some of these d = 1 symmetries may be 
related to (appear as products of) the local symmetries 
discussed next, depending on the topology enforced by 
boundary conditions. Some others are fundamental and 
cannot be expressed in terms of those local symmetries. 

d = symmetries: 
For the models under consideration, local, also called 
gauge, d = symmetries are associated with the ele- 
mentary loops (or plaquettes) P of the wires, see Fig. 1 
for an example. That is, when considering the supercon- 
ductors as point nodes, the links I form a network with 



minimal closed loops P. The associated local symmetries 
are given by 



Gi 



n( 

leP 



iciicu), Gp = t. 



(14) 



Repeating the considerations of (a) and (b) above, we 
see that, for any elementary plaquette P, the product of 
Majorana Fermi operators in Eq. (14) commutes with 
i^Mj since it shares an even number (possibly zero) of 
Majorana fermions with any bond in the Hamiltonian. 
By multiplying operators Gp for a collection of plaque- 
ttes P that, together, tile a region bounded by the loop 
r, it is readily seen that this product is also a symmetry, 
as in standard theories with gauge symmetries. 

The symmetries above lead to non-trivial conse- 
quences: 

(A) By virtue of Elitzur's theorem"^^ and its 
d > generalizations'"'"^^ all non-vanishing correlators 
(riaes ''") ■'^ith S a set of sites a must be invariant under 
all of the symmetries of Eqs. (13,14). That is, d = 0, 1- 
gauge-like symmetries cannot be spontaneously broken. 
As we alluded to earlier, one consequence of the non-local 
symmetries such as the d = 1 symmetries of Eq. (13) is 
the existence of TQO.24.37 

(B) Bounds on autocorrelation times. As a conse- 
quence of the d = 1 symmetries of Eq. (13), and the 
aforementioned generalization of Elitzur's theorem as it 
pertains to temporal correlators,-'® the Majorana Fermi 
system will exhibit finite autocorelation times regard- 
less of the system size. Of course, for various realiza- 
tions of dynamics and geometry of the disorder, differ- 
ent explicit forms of the autocorrelation times can be 
found. For instance, by use of bond-algebras, Kitaev's 
toric code model is identical to that of a classical square 
plaquette model as in Ref. 40. Similarly, Kitaev's toric 
code model"^ can be mapped onto two uncoupled one- 
dimensional Ising chains. '^'''^''^ Different realizations of 
the dynamics can lead to different explicit forms of r 
in both cases, however, finite autocorrelation times are 
found in all cases (as they must be). Similarly, more 
general than the exact bond algebraic mapping and di- 
mensional reductions that we find here, by virtue of d = 1 
symmetries of Eq. (13), autocorrelation functions involv- 
ing Majorana fermions on a line I must be bounded by 
corresponding ones in a d = 1 dimensional system.'® 



IV. DUALITIES AND SPIN REALIZATIONS OF 
ARBITRARY MAJORANA ARCHITECTURES 

In this Section, we provide two spin duals to the inter- 
acting Majorana system described by the effective Hamil- 
tonian Hm of Eq. (5) on arbitrary lattices/networks. 
This applies to finite or infinite systems and for arbi- 
trary boundary conditions. These two dual systems are 
(1) quantum Ising gauge theories for D = 2 systems, and 
more general spin gauge theories in higher-dimensions, 



and (2) a family of transverse-field Ising models with an- 
nealed disorder in the exchange couplings (each model 
representing a single gauge sector of -ffivi)- The dualities 
will be established in the framework of the theory of bond 
algebras of interactions, '^^''^^ as it applies to the study of 
general dualities between many-body Hamiltonians. The 
general bond algebraic method relies on a comparison of 
the algebras, in the respective two dual model, that arc 
generated by the corresponding local interaction terms 
(or bonds) in these theories. ^^'^^"■^'^ For the problem at 
hand, the Hamiltonian Hm is built as the sum of two 
sets of Hcrmitian bonds 



JQ1Q2, Vr, 



(15) 



where I and r are links and sites of the network sup- 
porting Hm (Vr was defined in Eq. (6)). In this paper, 
wc will only consider the bond algebra Am generated by 
these bonds. Wc can then obtain dual representations 
of Hm by looking for alternative local representations of 
Am- But first we have to characterize Am in terms of 
relations. 

The problem of characterizing a bond algebra of inter- 
actions is simplified by several features brought about by 
physical considerations of locality. The first consequence 
of locality is that interactions arc sparse^ meaning that 
each bond in any local Hamiltonian commutes with most 
other bonds and is involved in only a small number of 
relations (or constraints) that link individual bonds to 
one another. Hence the number of non-trivial relations 
per bond is small. The second consequence is that re- 
lations in a bond algebra can be classified into intensive 
and extensive, and most relations are intensive. We call 
a relation intensive if the number of bonds it involves is 
independent of the size of the system, and extensive if the 
numbers of bonds it involves scales with the size of the 
system. Since extensive relations could potentially lead 
to unphysical non-local behavior, they are typically few 
in number and may reflect the topology of the system 
regulated by the boundary conditions, as we will illus- 
trate repeatedly in this paper. As there are {2zr) Ma- 
jorana modes (or, cquivalcntly, Zr fermionic modes) per 
grain, the Majorana theory of Eq. (5) and the algebraic 
relations listed above are defined on a Hilbert space of 
dimension dimT^M = 2^'^^^ 

Next, we characterize the bond algebra Am as the first 
step toward the construction of its spin duals. The in- 
tensive relations are: 



1. for any r and I 

2. for r, r' € I, 

{Vr^iciicn} = = {Vr'^iciicn}, 

3. for r e ii, i = 1, 2, • • • , g^, 

{7'r,«Q,iQ,2} == 0. 



(16) 



(17) 



(18) 



Thus in the bulk, or everywhere for periodic boundary 
conditions, each island anticommutes with four qr (the 
coordination of r) links, and each link anticommutes with 
two islands. The presence or absence of extensive re- 
lations depends on the boundary conditions. For peri- 
odic (toroidal) or other closed boundary conditions (e.g., 
spherical), we have one extensive relation 



JjT'r =aJ|(JcnQ2), a = ±1, 



(19) 



since each Majorana Fermion operator appears exactly 
once both on the left and right-hand side of this equa- 
tion, but not necessarily in the same order. The constant 
a adjusts for the potentially different orderings, and the 
overall powers of i on each side of the equation. Notice 
that Y\^ Vr — Q 'iS the global Z2 symmetry operator. In 
contrast, for open or semi-open (e.g., cylindrical) bound- 
ary conditions, the islands on the free boundary have 
Majorana fermions that are not matched by links (that 
is, that do not interact with Majoranas on other islands). 
Hence the product 



wvA n(*^"^'2) =B 



(20) 



reduces to the product B of these Majoranas on the free 
boundary. The operator B may or may not commute 
with the Hamiltonian, depending on the details of the 
architecture at the boundary, sec Fig. 4, but cither way 
Eq. (20) does not represent an extensive relation in the 
bond algebra (rather it just states how to write a partic- 
ular operator as a product of bonds). If [Hm,B] = 0,-8 
represents a Z2 boundary symmetry independent of the 
local symmetries. 





FIG. 4. Two architectures with open boundary conditions. 
In either case, the operator B of Eq. (20) is the product of 
all the uncoupled Majoranas on the boundary indicated by 
open circles, but [Hu,B] — only for the system shown in 
the panel on the left. 



A. Duality to quantum Ising gauge theories 

In this Section we describe a duality relating the 
Hamiltonian H\j\ to a system of S = 1/2 spins. The 
spin degrees of freedom are placed on the (center of the) 
links of a network identical to the one associated to Hu , 
and are described by Pauli matrices erf, of , af . The goal 
is to introduce interactions among these spins that sat- 
isfy the same algebraic relations as the bonds of Hu ■ Let 
us introduce the Hermitian spin bond 



{l\rel} 



(21) 



For example, for the special case of the square lattice 
discussed in the Introduction, 

"Pr = <^tyi,^K. r^hMMM- (22) 

The set of spin bonds 

af, Vr, (23) 

satisfy the following intensive relations 

1. for any r and I 

{aff = 1 = {Vrr, (24) 

2. for r, r' £ I, 

{Vr,cTf}^0^{V,,,an, (25) 

3. ioT r e li, I = 1, 2, • • • , Qr, 

{P.,af}=0, (26) 

everywhere for closed boundary conditions, and every- 
where in the bulk for open or semi-open boundary con- 
ditions. These relations are identical to the intensive re- 
lations for the bonds of Hm- In the Ising gauge theory, 
the bond algebraic relations listed above are defined on 
a space of size 2^''^^ . (That this is so can be easily seen 
by noting that there are Ni = ZrNr links each endowed 
with a spin 5 = 1/2 degree of freedom af.) As it so hap- 
pens, this Hilbert space dimension is identical to that of 
the Majorana system of Hm- Putting all of the pieces 
together, we see that the spin Hamiltonian 



-f^QIG = - 2^ JlC^i - / ^^r 



(27) 



is unitarily equivalent to Hm , provided the extensive rela- 
tions are matched as well. For open or semi-open bound- 
ary conditions, the same follows provided that the inten- 
sive relations on the boundary also properly match. In 
the following, we focus on periodic boundary conditions 
(of theoretical interest in connection to TQO), and leave 
the discussion of open boundary conditions (of interest 



for potential experimental realizations of these systems) 
to the Appendix A. 

As just explained, the mapping of bonds 



»QiQ2 i-> o-i 



Vr ^Vr 



(28) 



preserves the intensive algebraic relations. In particular 
it maps the local symmetries of Eq. (14) to local sym- 
metries of Hq\q, 



Gp = ]J(icnQ2) '-^11'^'' = ^•S'.-P- 



(29) 



ieP 



leP 



To assess the effect it has on the extensive relation of Eq. 
(19) (and the global symmetry), notice that (for periodic 
boundary conditions) 



n(*cnQ2) '^11'^'' = *3'S 



(30) 



with [Qst Hq\q] a global symmetry of Hq]q, and 

l[rr^Y[rr = t. (31) 

r r 

It follows that, as it stands, the mapping of bonds of 
Eq. (28) is a correspondence, but not an isomorphism of 
bond algebras. The simplest way to convert it into an 
isomorphism is to modify one and only one of the bonds 
Vr of the spin model at some arbitrary site rg, so that 



n < 

{l\rael} 



(32) 



[a was defined in Eq. (19)) while for any other site 
r 7^ To, Vr remains unchanged. The introduction of 
this modified bond does not change the intensive rela- 
tions, since Qs commutes with every bond (original or 
modified). Moreover, 






(33) 



and the extensive relation of Eq. (19) is now, with the 
modified definition oiVra-, preserved, since (p? = 1) 



Xl^'t^^Vr, X{Vr 



(34) 



Hence there is a unitary transformation lA^ such that 

U^HuUl = i/QiG, (35) 

with Hq\q containing the single modified bond Vro ■ 

In the duality between the systems of Eqs. (5, 27), 
the dimensions of the their Hilbert spaces are identical. 
Since we count two Majorana modes (or, equivalently, 
one Fermionic mode) per link, the Hamiltonian Hu is 
defined on a Hilbert space of dimension dim^M = 2^', 



with Ni denoting the total number of hnks in the net- 
work. On the other hand, the spin system has one spin 
S = 1/2 degree of freedom per hnk, hence the dimension 
of the Hilbert space on which Hq\q is defined is also 2^' . 
Notice that the need to introduce the modified bond Vrg 
in the dual spin theory is irrelevant from the point of ex- 
ploiting the duality to study the ground-state properties 
of Hm (or viceversa, to study the ground-state proper- 
ties of Hq]q), since for finite systems the ground state 
\fi) must satisfy Qs\^) = |^)- The ease with which we 
established the duality between Majorana systems and 
QIG systems for general lattices and networks illustrates 
how efficient the bond-algebraic construct is. 



The duality just described is extremely general, valid 
in particular for any number of space dimensions D. In 
the following we describe explicitly one particularly im- 
portant special instance, that of Z) = 2. On a square 
lattice, the Hamiltonian Hm simplifies to 



Hh 



\Jl{iciiCi2) - y^/lrQiiC;3iC;32Q4 



(36) 



where li, 12,13,14 are shown in Fig. 1. This generalizes 
the Hamiltonian considered in Ref. 19 only in that inho- 
mogeneous couplings are allowed. The dual spin (finite- 
size) system is (ro € io,i, '0,2^0,3^0,4) 



r^ro I 



(37) 



r 



that we recognize as the standard, D = 2, Z2 Ising gauge 
theory, •^^ up to the modified bond at Tq {Qs = Yii '^i ^^^^ 
a is determined according to Eq. (19)), see Fig. 5. 




FIG. 5. Duality to a D = 2 Z2 quantum Ising gauge theory, 
where spins are represented as crosses. Hamiltonian Hm of 
Eq. (36) represents a particular fermionization of the Hq\c, 
gauge theory. 

Hence, we may regard the Hamiltonian of Eq. (36) as 
an exact fermionization of the Z2 Ising gauge theory with 
periodic boundary conditions (and one modified bond). 
It is interesting to compare this fermionization with a 
slightly different one^^ that exploits the Jordan- Wigner 
transformation in the limit of infinite size. This approach 
yields the Majorana Hamiltonian^^ (in our notation) 

Hfss = -^Jiiiciicu) - YKci^2Cir,iCi^2Ci2i, (38) 



where Ii,l2,l5,le are shown in Fig. 1. The two-body in- 
teraction Cij2Q5iQ62Q2i *'5 different than the two-body in- 
teraction in Hm, since it involves three different islands, 
see Fig. 6. Hence, disregarding boundary conditions, 
we see that the I12 quantum Ising gauge theory admits 
rather different but equivalent fermionizations. As ex- 
pected, the bonds in -ffpss satisfy intensive relations iden- 
tical to those already discussed for Hm and Hq\q. 




FIG. 6. Jordan- Wigner fermionization of the Z2 Ising gauge 
theory realizes a theory of Majorana fermions H^ss, with two- 
body interactions between Majoranas (shown as trapezoids) 
on three different islands. Notice that, unlike the intra-island 
two-body interactions of Hm, two neighboring two-body in- 
teractions -ffpss share a Majorana operator. 

Thus far, we focused on periodic boundary conditions. 
We now remark on other boundary conditions. When 



antipcriodic boundary conditions are imposed in a net- 
work with an an outer perimeter that includes twice an 
odd number of hnks, the right-hand side of Eq. (31) is 
replaced by —1. The union of both cases (periodic and 
antiperiodic) for a system having a twice odd perimeter 
spans all possible values of the product Ylr T^r- Thus, for 
these systems in the case of periodic boundary conditions, 
the spectrum of the Majorana system can be mapped to 
the union of levels found for the QIG systems for both 
periodic or antiperiodic boundary conditions. In terms 
of the corresponding partition functions, we have that 



= Z, 



QIG. periodic 



+ z, 



QIG, antiperiodic- 



(39) 



B. 



Duality to annealed transverse-field Ising 
models 



We next derive, in a similar spirit, a duality between 
the general architecture Majorana system Hm and an- 
nealed transverse-field Ising models. The number of an- 
nealed disorder variables in these systems (along with the 
number of sites TV,.) determines the size of the Hilbert 
space on which the Ising models are defined. With an 
eye towards things to come, we note (as we will re-iterate 
later on) that the duality that we will derive in this Sec- 
tion will furnish an example in which the Hilbert space 
dimensions of two dual systems need not be identical to 
one another. Generally, dualities are unitary transforma- 
tions between two theories up to trivial gauge redundan- 
cies that do not preserve the Hilbert space dimension. '^^ 
That is, dualities are isometrics. 

To define the annealed transverse field Ising systems, 
we place an S = 1/2 spin on each site r, cr^ , cr|f , cr^ , of 
the network associated to Hu, and a classical annealed 
disorder variable "qi = ±1 on each link I. Then we can 
introduce the set of Hermitian spin bonds 



rjia^a^ 



r,r' e I. 



(40) 



If we specialize to periodic boundary conditions, these 
bonds satisfy a set of intensive relations identical to the 
ones discussed in the two previous Sections, together with 
one new relation absent before and listed last below: 



1. for any r and I 



(r/id^cr^,) , 



2. for r,r' £ I 

{cr^,riia^a^,} = = {<,, Tyjcr^cr^,}, 

3. for r G ii, i = 1, 2, • • • , qr, 

4. for any elementary loop P in the network, 

leP leP 



(41) 



(42) 



(43) 



(44) 



The constraint of Eq. (44) holds true for any closed loop. 
For this reason, and others related to TQO, it is impor- 
tant to clarify the meaning of elementary loop. 

Loops in the network that share some links can be 
joined along those links to obtain another loop or sum of 
disjoint loops. This means that the set of all loops has a 
minimal set of generators from which we can obtain any 
loop or systems of loops by the joining operation just 
described. We call the loops in an arbitrary but fixed 
minimal generating set elementary loops. In this way, we 
obtain a minimal description of the constraints embod- 
ied in Eq. (44). It is not obvious a priori whether one 
should classify these constraints (that is, relations) as in- 
tensive or extensive. This depends on the topology of the 
system. If the system is simply connected, every loop is 
contractible to some trivial minimal (that is, of minimal 
length) loop, and hence we can choose minimal loops as 
elementary loops. These loops afford an intensive char- 
acterization of the constraints embodied in Eq. (44). If 
on the other hand the system is not simply connected, 
as for periodic boundary conditions, the generating set 
of elementary loops will include non-contractible loops, 
and the length of some of these non-contractible loops 
may scale with the size of the system. Consider, for ex- 
ample, the spin bonds of Eq. (40) on a planar network 
on the torus and on a punctured infinite plane. Both 
networks fail to be simply connected, but only the torus 
forces some of the constraint of Eq. (44) to be extensive, 
because its two non-contractible loops must scale with 
the size of the system. 

For periodic boundary conditions, there is one exten- 
sive relation satisfied by the bonds of Eq. (40), 



n('/'^rf^r') ='/!: 



with 



U^^ 



±1, 



(45) 



(46) 



that may or may not be independent of the relations of 
Eq. (44), depending on the details of the network. In 
the following, we will treat it as an independent relation, 
since it does not affect our results if it turns out to be 
dependent. 

It follows that the mapping of bond algebras 



iciicu H' rjialal, 



Vr-^al 



(47) 



preserves every local anticommutation relation. Hence 
the Hamiltonian theory 

H^^{m} = -Y, Mvi<J>r') - Y. ^-""r. (48) 



obtained from applying this mapping to H^ will be 
shown to be dual to Hu. (see Fig. 7). The Hilbert space 
on which the theory of Eq. (48) is defined is of size 
dim^Ai = 2^^"*"^'' where N^. is the number of supercon- 
ducting grains and N^, the total number of rji fields. 
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FIG. 7. Duality to an annealed transverse-field Ising model, 
in the particular D = 2 case. Spins S = 1/2 are located at 
the vertices r of the square lattice and classical Z 2 fields rji 
at the links I (indicated by a dash). 



The proposed duality raises an immediate question: 
What are the features of Hm that determine or at least 
constrain the classical fields rji? As we will see, the an- 
swer lies in the local and gauge-like symmetries that Hm 
possesses and Hm lacks. To understand this better, we 
need to study the effect this mapping has on relations 
beyond local anticommutation. Let us consider first its 
effect on the extensive relation of Eq. (19). We have that 

Q = l[Vr^l[<7-r^Qs, (49) 

r r 

a JJ(iQiQ2) i-^ a]^(?/icr^cr^') = aijt. (50) 

{ I 

As for periodic boundary conditions the left-hand sides of 
Eqs. (49) and (50) represent the same operator, but the 
right-hand sides arc different operators, the mapping as 
it stands docs not preserve the relation of Eq. (19). We 
know of a solution to this shortcoming from the previous 
Section. If we modify one and only one bond placed on 
some fixed but arbitrary link Iq to read 



mp]i,,a^(j^,Qs, 



(51) 



then 



a]^(iQiQ2) i-> 
I 

a^(TOo<<')Qs n ('?'^>r') = Qs, (52) 

as required by Eq. (19). 

The presence of the modified bond at Zq introduces a 
new feature into the discussion leading to Eq. (44). Now 
we have that, for any elementary loop P, 

i ^Ui^pVi, if lo^P, 

HiVK^'rO^l (53) 

leP { ctrjQs HieP 'ni if '0 e P- 



If we consider the role of the elementary loops P in the 
Majorana system Hy\^ and consider the mapping of Eq. 
(40), we see that the local symmetries (see Section III) 



Gp = W_{iciici2) 



(54) 



JeP 



of Hu are mapped to one of the two possibilities listed in 
Eq. (53), showing that as it stands, the mapping of Eq. 
(40) is still not an isomorphism of bond algebras. The 
problem is that a large number of distinct symmetries are 
being mapped either to a trivial symmetry (a multiple of 
the identity operator), or a multiple of the global Z2 
symmetry Qs of the annealed Ising model. We can fix 
this problem by decomposing the Hamiltonians Hy\ and 
Hi\\ into their symmetry sectors, where the obstruction 
to the duality mapping disappears. Thus we are able to 
establish emergent dualities, '^^'^^ that is, dualities that 
emerge between sectors of the two theories. 

The sector decomposition is simple for H/\\, that has 
only one symmetry Qs^ with eigenvalues qs = ±1. Then 
we can decompose the Hilbert space T^ai as 



qs=±l 



(55) 



so that if Aqg is the orthogonal projector onto T-Lq^ , then 



A„, = 5b A, 



iSi^qs 



qs- 



(56) 



For Hm, since its symmetries form a commuting set, 
one can simultaneously diagonalize them and break the 
Hilbert space Hm into sectors labelled by the symmetries' 
simultaneous eigenvalues, q = ±1 for the global symme- 
try and Tp = ±1 for the loop symmetries: 



"H-M — (37 ^(;,{rj=}- 
9,{rp} 



(57) 



The Hamiltonian Hm is block-diagonal relative to this 
decomposition, and, if A^ {pp} is the orthogonal projector 
onto the subspace Hg^{rp}, we have that 



Q^q,{rp} =<l^q,{rp}, 
GpAq^^rp} = TpAq^irp}, 



(58) 
(59) 



for any elementary loop P. 

The problem now is to decide which choice of sectors 
will make the projected Hamiltonians HMAq^Yp} and 
H/wAq^ dual to each other. From Eqs. (49), and (53), 
we obtain the relations 

q^qs, (60) 

f Uiepli^ if 'o^-P, 

Tp = ^ (61) 

[ ams rijep Vi if ^0 e P, 

which allow us to connect the two theories 

UdHuAq.iFpM = HM{vi}^qs, (62) 
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where the unitary transformation C/d implements an 
emergent duality that holds only on the indicated sectors 
of the two theories. 

The dual spin representation of Hm projected onto the 
gauge-invariant sector q = l,{rp = 1} is given by the 
inhomogeneous Ising model {rji = 1 on every link) 



H, 



Al 



{1} = - XI JlKK' - Yl ^^^r^ 



(63) 



and is known as a gauge-reducing duality.^^ For the spe- 
cial case of the square lattice and homogeneous couplings, 
one would expect that this sector contains the ground 
state of Hm . This latter result was derived, using meth- 
ods very different to ours, in Ref. 19. 



C. Physical Consequences 

We have by now seen, on general networks in an ar- 
bitrary number of dimensions, that ordinary quantum 
Ising gauge theories (and their generalizations) and an- 
nealed transverse-field Ising models arise from the very 
same Majorana system when it is dualized in different 
ways. Therefore, by transitivity. 



Hq\g < — > H/\\. 



(64) 

This correspondence leads to several consequences. In its 
simplest incarnation, that for D — 2 Majorana networks, 
this duality connects, via an imaginary-time transfer 
matrix (or T-continuum limit) approach, '^^'^^ disordered 
D = 3 classical Ising models to D = 3 classical Ising 
gauge theories. In its truly most elementary rendition 
among these planar networks, that of the square lattice, 
the duality of Eq. (64) implies that the effect of the bi- 
modal annealed disordering fields t;; = ±1 is immaterial 
in determining the universality class of the system. This 
is so as the standard random transverse-field Ising model 
on the square lattice 



Hrtpm = — 2_^ Jl<yr'^r' ^ /_^ hr<7^ 



(65) 



(i.e., Eq. (48) in the absence of annealed bimodal dis- 
order) similarly maps, via a transfer matrix approach, 
onto a corresponding classical Ising model on a cubic lat- 
tice. The uniform transverse-field Ising model (that with 
uniform Ji and hr) maps onto the uniform D = 3 Ising 
model. Thus, in this latter case, the extremely disor- 
dered system with annealed random exchange constants 
exhibits the standard D = 3 Ising type behavior of uni- 
form systems. 

By the dualities of Sections IV A and IV B, general 
multi-particle, or multi-spin, spatio-temporal correlation 
functions in different systems can be related to one an- 
other. In particular, by Eq. (28) relating the Majorana 
system with the quantum Ising gauge theory, the two 
correlators 

(n^^W(*QiQ2)(t')> = {HVritWnt')) (66) 



are equal. Thus, if certain correlators (e.g., standard 
static two-point correlation functions, autocorrelation 
functions, or four-point correlators such as those preva- 
lent in the study of glassy systems)''^ appear in the spin 
systems, then dual correlators appear in the interacting 
Majorana system with identical behavior. An exact du- 
ality preserves the equations of motion, and so the dy- 
namics of dual operators are the same."^^ Similarly, by 
the duality of Eq. (35), the phase diagrams describing 
the Majorana networks arc identical to those of quantum 
Ising gauge systems. In instances in which the quantum 
Ising gauge theories have been investigated, the phase 
boundaries in the Majorana system may thus be mapped 
out without further ado. 

Lattice gauge theories with homogeneous cou- 
plings, i.e., uniform lattices, have been investigated 
extensively. ■^^'^"'^ As we alluded to above, it is well appre- 
ciated that the quantum Ising gauge theory on a square 
lattice can be related, via a Feynman mapping, to an 
Ising gauge theory on the cubic lattice with the classical 
action 



S\c 



-kJ2vp. 



(67) 



rJ 



r,l 



The latter has a transition"'^ a.t K ~ K^ — 0.761423, a 
value dual^^ to the critical coupling (or inverse critical 
temperature when the exchange constant is set to unity) 
of the D — 3 classical Ising model with nearest neighbor 
coupling, Kc = 0.2216595. Similar transitions between 
a confined (small K) to a deconfined (large K) phases 
appear in general uniform coupling lattice gauge theories 
with other geometries. Phase transitions mark singulari- 
ties of the free energy, that are always identical in any two 
dual models. ^^ In our case of interest here, by the corre- 
spondence of Eq. (35), identical transition points must 
thus appear in the dual Majorana theories. In particular, 
the transition points in the Majorana system are imme- 
diately determined by their dual spin counterpart. More 
precisely, the Majorana uniform network depicted in Fig. 
1 displays a quantum critical point of the D ~ 3 Ising uni- 
versality class at {J/h)c = — 2_ftrc/lntanh A'c = 0.29112. 
In theories with sufficient disorder (e.g., quenched ex- 
change couplings, fields, or spatially varying coordina- 
tion number), rich behavior such as that exemplified by 
spin glass transitions or Griffiths singularities^^ may ap- 
pear. According to Eq. (35), in architectures with non- 
equidistant superconducting grains of random sizes, the 
effective couplings {Ji} and {hr} are not uniform and 
may lead to spin glass, Griffiths, or other behavior when- 
ever the corresponding dual gauge theory exhibits these 
as well. We note that the random transverse field Ising 
model of Eq. (65) is well known to exhibit a (quantum) 
spin glass behavior. ^^'^^ If and when it occurs, glassy (or 
spin-glass) dynamics in the annealed or gauge spin sys- 
tems will, by our mapping, imply corresponding glassy 
(or spin-glass) dynamics in the Majorana system as well 
as interacting electronic systems (leading to electron glass 
behavior) . The disordered quantum Ising model was em- 
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ployed in the study of the insulator to superconducting 
phase transition in granular superconductors.'*^ Numer- 
ous electronic systems are indeed non-uniform^^ and/or 
disordered. ^° 



SPIN DUALS TO SQUARE LATTICE 
MAJORANA SYSTEMS 



Thus far, we provided a systematic analysis of symme- 
tries and dualities for Majorana systems supported on 
networks in any number of spatial dimensions. It is in- 
structive to consider particularly simple architectures as 
these highlight salient features and, on their own merit, 
provide new connections among well studied theories. In 
what follows, we will focus on the square lattice super- 
conducting grain array of Fig. 1, and some honeycomb 
and checkerboard lattice spin dual models. 



A. The XXZ Honeycomb Compass Model 

The Majorana system Hm of Eq. (5) in a square lat- 
tice is dual to a very interesting spin Hamiltonian on the 
honeycomb lattice, see Fig. 8. The dual spin model may 
be viewed as an intermediate between the classical Ising 
model on the honeycomb lattice (involving products of a 
single spin component (tr^) between nearest neighbors) 
and Kitaev's honeycomb model, *^ for which the bonds 
along the three different directions in the lattice are re- 
spectively pairwise products of the three different spin 
components. This particular spin Hamiltonian, which 
we dub XXZ honeycomb compass modeX is described by 



H, 



XXZh — 






J I crier 



r+ei 



-J2k 



r^e. 



.(68) 



where each S = 1/2 is located on the vertices r of a 
honeycomb lattice, and ap^ are the corresponding Pauli 
matrices. The qualifier "non-vertical links" alludes to 
the two diagonally oriented directions of the honeycomb 
lattice while "vertical links" are, as their name suggests, 
the links parallel to the vertical direction in Fig. 8. The 
unit vector ei points along the diagonal link I and may be 
oriented along any of the two diagonal directions. The 
XXZ honeycomb compass model exhibits local symme- 
tries associated with every lattice site r. 



G: 



XXZh 



(69) 



Similarly, the XXZ system exhibits d = 1 symmetries of 
the form 



iXXZh 






(70) 



associated with every non- vertical contour £ (i.e., that 
composed of the diagonal non- vertical links) that circum- 
scribes one of the toric cycles. 

We provide, in the left-hand panel of Fig. 8, a sim- 
ple schematic of the topology of the honeycomb lattice - 




FIG. 8. The brick-wall planar orbital compass model 
(shown on the left) can be seen as a simpler relative of the 
XXZ honeycomb compass model, by placing it on a honey- 
comb lattice as shown on the right. 



that of a "brick- waU lattice" .^°'^^ The brick- waU lattice 
also captures the connections in the honeycomb lattice. 
It is formed by the union of the highlighted vertical (red) 
and horizontal (green) links in the left-hand side Fig. 8. 
The brick-wall lattice can be obtained by "squashing" 
the honeycomb lattice to flatten its diagonal links while 
leaving its topology unchanged in the process. In the 
brick- wall lattice, ei simply becomes a unit vector along 
the horizontal direction. As can be seen by examining 
either of the panels of Fig. 8, the centers of the vertical 
links of the honeycomb (or brick- wall) lattice form, up to 
innocuous dilation factors, a square lattice. As is further 
evident on inspecting Fig. 8, between any pair of centers 
of neighboring vertical (red) links, there lies a center of 
a non-diagonal (green) link. This topological connection 
underlies the duality between the Majorana model on 
the square lattice and the XXZ honeycomb compass spin 
model. We explicitly classify the bonds in the Hamilto- 
nian of Eq. (68) related to the two types of geometric 
objects: 

1. Bonds of type (i) are associated with the products 
Wr'^r+e } "^^ diagonal links of the lattice. They 
each anticommute with two 

2. Bonds of type (ii), affiliated with products 
{'^r'^r+e } '^^ ^'^^ vertical links. Each one of these 
bonds anticommute with four bonds of type (i). 

We merely note that replacing the bonds of the Ma- 
jorana model on a square lattice, as they appear in the 
bond algebraic relations (1-3) of Section IV, by the ones 
above leads to three equivalent relations that completely 
specify the bond algebra of the system of Eq. (68). As 
we have earlier seen also the quantum Ising gauge theory 
of Eq. (27) and the annealed transverse-field Ising model 
of Eq. (48) have bonds that share the same three ba- 
sic bond algebraic relations. Thus we conclude that the 
XXZ honeycomb compass model is exactly dual to the 
quantum Ising gauge theory of Eq. (27) on the square 
lattice. In its uniform rendition (with all couplings Ji 
and fields hr being spatially uniform) the XXZ honey- 
comb compass system lies in the 3D Ising universality 
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class. Similarly, many other properties of the XXZ hon- 
eycomb compass model can be inferred from the heavily 
investigated quantum Ising gauge theory. 

The duality between the XXZ honeycomb compass 
model and its Majorana system equal on the square lat- 
tice affords an example of a duality in which the Hilbert 
space size is preserved as we now elaborate. The XXZ 
theory of Eq. (68) is defined on a Hilbert space of 
size dimHxxzh ~ 2^''' where A^hi is the number of sites 
on the honeycomb lattice while that of the Majorana 
model of Eq. (5) was on a Hilbert space of dimension 
dimT^M = 4^''. Now, for a given number Nr of vertical 
links on the honeycomb lattice, we have the same number 
of bonds of type (i) and (ii) as we had in the Majorana 
system while having A^hi = 2Nr lattice sites. 



B. Checkerboard model of (p 
grains 



ip) superconducting 



In Ref. 51, Xu and Moore, motivated by an earlier 
work of Moore and Lee,^^ proposed the following spin 
Hamiltonian 



HxM — — /_^{hr '^r 






(71) 



to describe the time-reversal symmetry breaking charac- 
teristics in a matrix of unconventional p-wave granular 
superconductors on a square lattice. In writing Eq. (71), 
we employ a shorthand 



Daf = 



hei+ea^ 



(72) 



to denote the square lattice plaquette product, where ei 
and 62 denote unit vectors along the principal lattice 
directions. It is important to emphasize that the spins 
a^'^ in Eqs. (71), (72) are situated at the vertices r of 
the square lattice (not on the links (or link centers) as 
in gauge theories). The eigenvalues a^ = ±1 describe 
whether the superconducting grain located at the vertex 
of the square lattice r has a, {p + ip) or a (p — ip) order 
parameter. 

We show next that a D ^ 2 checkerboard rendition of 
the XM model which we denote by CXM (see Fig. 9) is 
dual to the Majorana system on the square lattice (which 
is, as we showed, dual to the XXZ honeycomb compass 
model and all of the other models that we discussed ear- 
lier in this work). This system is defined by the following 
Hamiltonian 



H, 



CXM 






E 






(73) 



a^l+^2— odd 



In this system, the plaquette operators Dct^ (with r = 
xiBi + X2e-2) appear in every other plaquette (hence 
the name "checkerboard"). These plaquettcs arc present 
only if xi + X2 is an odd integer as emphasized in Eq. 
(73). The model has the following local symmetries 



Gp^n 



(74) 



^(r 



62 



^ ( : 

n 



D 

^i : 



n 




n 



T-eP 



► ei 

FIG. 9. The checkerboard Xu-Moore (CXM) model of Eq. 
(73). The symmetry plaquettes P constitute half of all the 
plaquettes of the lattice, while the interaction plaquettes Dcr^ 
represent the other half. 



where P arc those plaquettcs appearing whenever a;i -|-a;2 
is an even integer. 

The proof of our assertion above concerning the 
duality of this system to the Majorana system of Eq. (5) 
when implemented on the square lattice is straightfor- 
ward and will mirror, once again, all of our earlier steps. 
We may view the Hamiltonian of Eq. (73) as coinprised 
of two basic types of bonds: 



1. Bonds of type (i) arc on-site operators {cr^} asso- 
ciated with local transverse fields. 

2. Bonds of type (ii) are the plaquette product op- 
erators {ncr^} of Eq. (72), for plaquettes whose 
bottom left-hand corner r is an "odd" site. 

The basic network structure underlying these bonds is 
simple and, apart from an interchange of names, identi- 
cal to that of the Majorana system on the square lattice 
of Fig. 1 as well as that of the XXZ honeycomb com- 
pass model of Fig. 8. To see this, we note that in the 
checkerboard of Fig. 9, the four-fold coordinated inter- 
action plaquettes generate, on their own, a square lattice 
grid. Between any two neighboring interaction plaque- 
ttcs on this square lattice array, there is a lattice site 
r (see Fig. 10). As in our earlier proof of the duality, 
we simply remark that replacing the bonds of the Ma- 
jorana model on a square lattice, as they appear in the 
bond algebraic relations (1-3) of Section IV, by the ones 
above leads to three equivalent relations that completely 
specify the bond algebra of the CXM system. The Ma- 
jorana and CXM models are thus dual to one another 
{Hu ■!->■ Hqxm) when their couplings are related via the 
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correspondence 



Ji ^ h: 



XM 



XM 



(75) 



Thus, the CXM model joins the fellowship of all other 
dual theories (with the same network connectivity) that 
wc discussed in this work (i.e., the Majorana, quantum 
Ising gauge, and annealed transverse field Ising models 
on the square lattice as well as the XXZ compass model 
on the honeycomb (or equivalent brick- wall) lattice). 




FIG. 10. The D = 2 checkerboard Xu-Moore (CXM) model 
is dual to the Majorana system in a square lattice as shown 
on the left. On the right, we rotate and redefine the lattice 
in a manner which highlights its connection to the quantum 
Ising gauge (QIG) theory of Eq. (27). 

On the right-hand half of Fig. 10, we pictorially illus- 
trate the connection between the CXM model and the 
quantum Ising gauge theory. The individual sites of the 
checkerboard lattice of Fig. 9 (the sites at which the lo- 
cal transverse fields are present) map onto links of the 
gauge theory (Section IV A). Similarly, the interaction 
plaquettes of the CXM model map into plaquettes of the 
quantum Ising gauge theory. Note, on the right, that as 
is geometrically well appreciated, the four center-points 
of the individual links on the square (gauge theory) lat- 
tice can either circumscribe interaction plaquettes of the 
gauge theory or may correspond to four links that share a 
common endpoint that do form a "star" configuration.'^^ 
In particular, by its duality to the quantum Ising theory, 
the CXM rigorously lies in the 3D Ising universality class 
when the couplings J^ and h/^^ are spatially uniform. 
For a given equal number of bonds in both the Majorana 
system and the CXM theory, it is readily seen that the 
Hilbert space dimensions of both theories are the same, 
dim^M = dimHcxM. 



VI. SIMULATING HUBBARD-LIKE MODELS 
WITH MAJORANA NETWORKS 

The Dirac, fermionic, annihilation and creation oper- 
ators, {dr} and {d\.} respectively, can be expressed as a 
linear combination of two Majorana fermion operators. 
For example, if we are interested in two-flavor Dirac op- 



erators a possible realization is (see Fig. 1) 



dri = -jE,{ci^i + ici^2), dL = —7^{ci^i " ^c^z), (76) 



V2^ 



1 

71^ 



where r € Zi, ^2, '3, li- 

A system of interacting Dirac fermions (e.g., electrons) 
on a general graph can be mapped onto that of twice the 
number of Majorana fermions on the same graph, and 
each Dirac fermion is to be replaced by two Majorana 
fermions following the substitution of Eq. (76). Thus, 
any granular system of the form of Eq. (5) in which each 
grain r has qr = 2zr neighbors, can be mapped onto a 
Dirac fermionic system on the same graph in which on 
each grain there are z^. Dirac fermions. There are many 
possible ways to pair up the Majorana fermions in the 
system of Eq. (5) to yield a corresponding system of 
Dirac fermions. Equation (76) represents just one possi- 
bility. Another possible way to generate (spinless) Dirac 
fermions is 

di = ^(Qi + *cj2), dj = -7^{cii - icn)- (77) 

All of the spin duals that we derived for Majorana 
fermion systems hold, mutatis mutandis, for these sys- 
tems of Dirac fermions on arbitrary graphs. In this 
sense, dualities afford an alternative, flexible approach 
to fermionization that does not rely on the Jordan- 
Wigner transformation.^^ Most importantly, one can use 
these mappings to simulate models of strongly interact- 
ing Dirac fermions, such as Hubbard-like models, on the 
experimentally realized Majorana networks. In other 
words, one can engineer quantum simulators out of these 
Josephson junction arrays. 

As a concrete example, we consider the square-lattice 
array of Fig. 1 and transform, on this lattice, the Majo- 
rana system of Eq. (5) into a two-flavor Hubbard model 
with compass-type pairing and hopping. Based on our 
analysis thus far we will illustrate that this variant of the 
2D Hubbard model is exactly dual to the 2D quantum 
Ising gauge theory and thus lies in the 3D Ising univer- 
sality class. Consider the mapping of Eq. (76). With 
UrcF = dl,^dra ("" =t5i)i ^ Hubbard type term with on- 
site repulsion Ur becomes 



UriUr^ - l){nrl - 1) = Ur(Vr - 1), 



(78) 



akin to the second term of Eq. (5) with hr <-> Ur (up to 
an irrelevant constant). In what follows we assume that 
the network array of Fig. 1 has unit lattice constant. 

The Majorana bilinear that couples, for instance, the 
bottom most corner of the grain that is directly above 
r (i.e., site r -\- €2) to the top-most site of grain r (with 
thus a link I that is vertical) becomes 



-iJlCl^lCl22 



Ji 



idl+dri)id, 



■t 

r+e2l 



He2i 



)• (79) 
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Similarly, for horizontal links I, the bilinear in the first 
term of Eq. (5) realizes pairing hopping terms involving 
only the ^ flavor of the fermions. Thus, the Hamilto- 
nian of Eq. (5) becomes a Hubbard type Hamiltonian 
with bilinear terms containing hopping and pairing terms 
between electrons of the up or down flavor for links I 
that are vertical or horizontal, respectively. Such a de- 
pendence of the interactions between the internal spin 
flavor on the relative orientation of the two interacting 
electrons in real-space bears a resemblance to "compass 
type" systems. ^^ Putting all our results together, the 
Dirac fermion Hamiltonian on the square lattice with pair 
terms of the form of Eq. (79) augmented by the on-site 
Hubbard type interaction term of Eq. (78) is dual to all of 
the other models that we considered thus far in this work. 
In particular, as such this interacting Dirac fermion (or 
electronic) system is not of the canonical non-interacting 
Fermi liquid form. Rather, this system lies in the 3D 
Ising universality class. 

The standard Hubbard model with SU{2) spin sym- 
metry, which up to chemical potential terms is given by 
(a = 1,2) 



^Hub — —t 2_^ idra'^r+e^cr + h.C.) 



+U^inrt-l){nri-l). 



(80) 



can be written as a sum of terms of the form of Eq. (78) 
augmenting many Majorana Fermi bilinear coupling sites 
on nearest neighbor grains (i.e., r and r ± e^). As we 
illustrate in Fig. 11, we label the four Majorana modes on 
each grain r as {cra}t=i- In terms of these, the Hubbard 
Hamiltonian becomes 



-H^Hub 



E 

a.a—l 



t y ^(CraCr+eaa+2 

r.a.a— 1.2 



Cr-feo, a^ra-f 2 j 



(81) 



Thus, the Hubbard Hamiltonian may be simulated via 
Majorana wires with multiple Josephson junctions. 

Appendix B describes the possible simulation of quan- 
tum spin S ~ 1/2 systems in terms of Majorana net- 
works. 



VII. CONCLUSIONS 

We conclude with a brief synopsis of our findings. This 
work focused on the interacting Majorana systems of 
Eq. (5) on general lattices and networks. By employing 
the standard representation of Dirac fermions as a linear 
combination of Majorana fermions, our results similarly 
hold for a general class of interacting Dirac fermion sys- 
tems on general graphs. Towards this end, we heavily 
invoked two principal tools: 

• The use of d-dimensional gauge-like symmetries 
that mandate dimensional reduction and TQO via 




FIG. 11. A labeling of the Majorana wire endpoints on the 
square lattice which we use here to explicitly represent the 
standard electronic Hubbard model in terms of Majorana op- 
erators. This is a different labeling than the one in Fig. 1. 



correlation function bounds.^* ^^ These symme- 
tries lead to bounds on the autocorrelation times. ^^ 



26,28-33 



as it, 



• The bond-algebraic theory of dualities 
in particular, pertains to very general dualities and 
fermionization^^''^^ to obtain multiple exact spin 
duals to these systems, in arbitrary dimensions and 
boundary conditions, and for finite or infinite sys- 
tems. 

Using this approach, we demonstrated that 

• The Majorana systems of Eq. (5), standard 
quantum Ising gauge theories (Eq. (27)) and, 
transverse-field Ising models with annealed bi- 
modal disorder (Eq. (48)) are all dual to one an- 
other on general lattices and networks. The duality 
afforded an interesting connection between heav- 
ily disordered annealed Ising systems and uniform 
Ising theories. The spin duals further enable us 
to suggest and predict various transitions as well 
as spin-glass type behavior in general interacting 
Majorana fermion (and Dirac fermion) systems. 
The representation of Dirac fermions via Majorana 
fermions enlarges the scope of our results. In par- 
ticular, as Eq. (78) makes evident, the standard 
on-site Hubbard term in electronic systems is ex- 
actly of the same form as that of the intra-grain 
coupling in the interacting Majorana systems that 
we investigated. We similarly represented the bi- 
linear in the Majorana model of Eq. (5) as a Dirac 
fermion form (Eq. (79)). Following our dualities. 
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on the square lattice, the interacting Dirac fcrmion 
(or electronic) Hamiltonian formed by the sum of 
all terms of the form of Eqs. (78, 79) is dual to the 
quantum Ising gauge theory and thus lies in the 3D 
Ising universality class, notably different from stan- 
dard non-interacting Fermi liquids; this non-trivial 
electronic system features Hubbard on-site repul- 
sion augmented by "compass" type hopping and 
pairing terms. We further showed how to simulate 
bona fide Hubbard type electronic Hamiltonians via 
Major ana wire networks. 

• Several new systems were introduced and investi- 
gated via the use of bond algebras: 

(1) the "XXZ honeycomb compass" model of Eq. 
(68) (a model intermediate between the classical 
Ising model on the honeycomb lattice and Kitaev's 
honeycomb model and, 

(2) a checkerboard version of the Xu- Moore model 
for superconducting (p + ip) arrays (Eq. (73)). 
By the use of dualities, we illustrated that both of 
these systems lie in the 3D Ising universality class. 

As evident in our work, all of the considerations nec- 
essary to attain these results were, to say the least, very 
simple by comparison to other approaches to duality that 
generally require far more involved calculations. In the 
appendices we discuss other connections between Majo- 
rana and spin systems. 
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Appendix A: Dualities in finite systems with open 
boundary conditions 

We have, so far, studied exact dualities for the Majo- 
rana system with the Hamiltonian Hm of Eq. (5) when 
subject to periodic boundary conditions. We focused on 
periodic boundary conditions these are pertinent to the 
theoretical study of TQO. In this appendix, we will con- 
sider exact dualities in the presence of open boundary 
conditions. In doing so, we will further study finite, even 
quite small, square lattices. It is useful to provide a pre- 
cise description of these finite dual spin systems as there 
is a definite possibility that this Majorana architecture 
may become realizable in the next few years. These dual- 
ities also allow us to illustrate the flexibility of the bond 
algebraic approach to dualities in handling a variety of 
boundary conditions exactly. As in the rest of this pa- 
per, the dualities we obtain are exact unitary equiva- 
lences. Thus, these dualities may be tested numerically 
by checking if the energy spectra of the two dual systems 
are indeed identical. 





FIG. 12. The spin dual of two superconducting islands. Each 
island maps to a plaquette interaction of the quantum Ising 
gauge theory, but such a mapping would not be compatible 
with matching dimensions of Hilbert spaces. Hence one of the 
lower plaquette is chopped to include only one spin. 



As illustrated in Section IV A, the effective Hamilto- 
nian Hm on the square lattice and in the bulk is dual to 
the ^2 lattice gauge theory. In this appendix, our task 
is to find the boundary terms that make the duality ex- 
act in the presence of open boundary conditions. Here 
we only consider dualities that preserve the dimension 
of the Hilbert space of the two theories. We thus fol- 
low two guiding principles: 1) in the bulk, the dual spin 
theory remains the Z2 lattice gauge theory, and 2) on 
the boundary, we introduce terms that preserve both the 
bond algebra and the dimension of the Hilbert space. Let 
us start with the simplest interacting case, that of two is- 
lands (grains) linked by one Josephson coupling, see Fig. 
12. In this case, the Hamiltonian of Eq. (5) reads 



Hm = -hciC2C3C4 — h'c^cecrcs — Jic^c^. 



(Al) 



This Hamiltonian acts on a Hilbert space of dimension 
dim'Hfvi = 2®/^ = 2*. Thus, the dual theory must contain 
four spins and some recognizable gauge interactions. The 
result is 



Hq\G = —h(7i — h CT]^ (T2 (T3 (74 — J(Ti . 



(A2) 



where the single spin af in the Hamiltonian stands for 
an incomplete plaquette. One can check that the bond 
algebra is preserved and the two spectra are identical. 

The next interesting case contains four superconduct- 
ing islands, see Fig. 13. In this case, dimHw = 2^^/^ = 
2®, and so the dual spin Hamiltonian, described diagram- 
matically in Fig. 13 contains eight spins, two complete 
and two incomplete gauge plaquettes. The situation be- 
comes more regular if we further increase the number 
of islands. For nine islands (dim^M = 2'^*'/^ = 2^*), the 
Majorana system maps to eighteen spins, three complete, 
and six incomplete plaquettes on the first and last row 
of the spin model. One can generalize this picture to 
L^ islands. Then the dual Z2 quantum Ising gauge the- 
ory will be represented by a scaled version of the right 
panel of Fig. 14, with 2i^ spins, and 2L incomplete pla- 
quettes (the product of only three spins a^ ) . The latter 
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FIG. 13. The spin dual for a configuration of four islands. 
The incomplete plaquettes represent two-spin interactions in 
the Hamiltonian. 
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FIG. 14. The spin dual for nine islands. Incomplete plaque- 
ttes represent three-spin interactions in the spin Hamiltonian, 
the product of the three spins a'^ closer to an incomplete green 
diamond. 



incomplete plaquettes are equal split between the top and 
bottom rows, i.e., L incomplete plaquettes are placed on 
the top row and L arc situated on the bottom row. 

Notice that there is no natural guiding principle to find 
the dual theory by a Jordan- Wigner mapping. The bond- 
algebraic method is the natural approach and can be tested 
numerically on finite lattices. 



Appendix B: Fermionization of S = 1/2 spin models 
in arbitrary dimensions 

Although not pertinent to our direct models of study 
(those of Eq. (5) and their exact duals), we briefly re- 
view and discuss, for the sake of completeness and general 
perspective, dualities of related quantum spin S = 1/2 
systems. General bilinear spin Hamiltonians can be ex- 
pressed as a quartic form in Majorana fermion opera- 
tors. The general nature of this mapping is well known 
and has been applied to other spin systems with several 
twists. Simply put, we can write each spin operator as 
a quadratic form in Majorana fcrmions. In the case of 
general two-component spin systems that we discuss now, 
the relevant Pauli algebra is given by the following on-site 
(r) constraints 



and trivial off-site (r 7^ r') relations, 



(B2) 



A dual Majorana form may be easily derived as follows. 
We consider a dual Majorana system in which at each 
lattice site r, there is a grain with three relevant Majo- 
rana modes. We label the three relevant Majorana modes 
(out of any larger number of modes on each grain) by 
{cr.,a}a=i- As Can bc readily seen by invoking Eq. (4), a 
representation that trivially preserves the algebraic rela- 
tions of Eqs. (Bl, B2) is given by 



a^ -i-> lCrlCr2, CTr "^ ICrlCr^. 



(B3) 



Equation (B3) is a variant of a well known mapping appli- 
cable to three component spins (as well as, trivially, spins 
with any smaller number of components). ^^'^"^ Equation 
(B3) may also be viewed as a two-component version of 
the mapping employed by Kitaev.^^ The Hilbert space 
spanned by an S" = 1/2 spin system on a lattice/network 
having N sites is dm\T-Lsp\n = 2^. By contrast, the 
Hilbert space of a general Majorana system with {m^} 
Majorana modes (nir > 3) at sites {r} is given by 
dim-HM = 2^-™'-/2. Thus, in this duality the Hilbert 
space is not preserved: each individual energy level of 
the spin system becomes 2^"^^ m^/2)-N ^qY(\ degenerate. 
Similarly, one-component systems (e.g., those involving 
only {cr^}) can be mapped onto a granular system with 
two Majorana modes per site. If there are two Majorana 
modes at each site r then such a mapping will preserve 
the Hilbert space size. 

For completeness, we now turn to specific spin systems 
related to those that we discussed in the main part of our 
article. In Section VB, we illustrated that the Majorana 
system of Eq. (5) (and all of its duals that we earlier dis- 
cussed in the text) can be mapped onto the Xu-Moore 
model^^ on the checkerboard lattice. Following our gen- 
eral discussion above, it is straightforward to provide a 
Majorana dual to the Xu-Moore model on the square lat- 
tice, Eq. (71). On the square lattice, the orbital compass 
model (OCM) and the Xu-Moore model of Eq. (71) are 
dual to one another. '^^'•^^^^^ We will assume the square 
lattice to define the xz plane. The anisotropic square 
lattice OCM^^'^^ is given by the Hamiltonian 



i?o 



CM 



(Jx-r<^r'^r 



J. 



z:r'^r'^r+B2 



). (B4) 



(a^ 



x\2 



(a^)2 = l, {a^,a^} = 0. 



(Bl) 



In Eq. (B4), we generalized the usual compass model 
Hamiltonian by allowing the couplings {Jx,z} to vary lo- 
cally with the location of the horizontal and vertical links 
of the square lattice (given by (r, r -I- 61^2) respectively). 
By plugging Eq. (B3) into Eq. (B4), we can rewrite this 
(as well as other general two-component spin bilinears) 
as a quartic form in the Majorana fermions. 

It may generally be feasible to use our formalism to 
simulate quantum spin models in terms of Majorana 
networks. Consider, for example, the simulation of a 
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FIG. 15. The transverse-field Ising model can be simulated 
by an architecture of nanowires with one wire per supercon- 
ducting island. 



transverse-field Ising chain 

7V-1 



N 



H\ = -J2 J-'^t'^t+i - Y. '^^''^ 



(B5) 



i=i 



with TV spins and open boundary conditions. In this 
case, it may be possible to use linear arrays with one 
nanowire per island to simulate this model and study, for 
instance, the dynamics of its quantum phase transition. 
The Hamiltonian H\ maps to the Majorana network 



N-l 



N 



Hk 



y^ JiCt,2Ci+i,i - i^^hid^i 



Ci.2, 



(B6) 



4 = 1 4=1 

after the following duality mapping 

<yt<yt+i H> iCi^2Ci+i,i , erf 1-)- iCi,iCi,2, (B7) 

see Fig. 15. 
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On each grain, the z^ non-intersecting nanowires link one 
half of the nanowire endpoints to the remaining half; 
there are {2zt.)\ / {2'^'^ Zr\) distinct ways for different pair- 
ings of the vertices. These non-intersecting nanowires can 
be placed in any way on a surface of the bulk supercon- 
ducting grain. For instance, in the square and triangular 
lattices, the regular arrangement of nanowires shown in 
Figs. 1, and 3 is only one among many others. 
It us useful at this point to recall some basic algebraic 
facts about Majorana fermions. Let us label the Majorana 
operators simply as Ci, i = 1, ■ • ■ , A^. Then, 



{ci,Cj} — 2(5i,. 



In general. 



(ci' 



■Cjv) 



= (-1) 



N(N-l)/2 



1, 



(B8) 



(B9) 



so that the eigenvalues of ci • • • C]v are ±1 or ±i, depending 
on A^. If A'^ is even, up to equivalence, there is only one 
irreducible representation (irrep) for these relations. It acts 
on a Hilbert space of dimension dimHM(A^ even) — 2^'^, 
and can be described in terms of Pauli matrices as 






•o-i). 



i=l,.- 
i=l-F 



JV 
2 ' ' 



N_ 
2 ' 



,A^. 



(BIO) 
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If N is odd, 



[Ci ■ ■ ■CN,Ci] = 0. 



(Bll) 



So in view of Eq. (B9), the irreps of Eq. (B8) are charac- 
terized by an extra constraint, 

C1---CN = al, (irrep for Af odd) (B12) 

with a = ±1, ±z depending on A'^ and the particular irrep. 
An exphcit irrep is afforded by 

c^ = (B13) 



o-i 



-.« 



• CTjV-l 



■••crrj 



JV+1 
2 ' 

A'. 



.A-- 



'(B14) 



Its dimension is d:YmHu{N odd) = 2''^"^'''^ 

^® A definite order is to be assigned on each grain. As will 
become evident, an odd permutation of the ordering of the 
Majorana operators in the product of Eq. (6), which leads 
by virtue of the Majorana algebra to a sign change, will 
not change the honA algebra (to be defined later) and thus 
none of our dualities. 

^^ In these and other general systems, not all d = 1 gauge-like 
symmetries are independent. If the Majorana system has 
a trivial homology (such as that of an infinite plane or a 
sphere), no d = 1 gauge- like symmetries appear: all sym- 
metries involving 0{L^) sites can be expressed in terms of 
a product of the local symmetries and thus are not fun- 
damental. By contrast, for a D — 2 dimensional system 
placed on a torus there are two closed toric cycles £ inde- 
pendent of the local symmetries. This explains the findings 
of Ref. 19 that the Majorana theory has TQO. Some man- 
ifestations of this phenomenon were noticed in numerical 
simulations of compact QED^^, though not recognized as 
a signal of the presence or absence of TQO. In numerous 
theories with plaquette and/or link interactions, the Euler- 
Lhuillier formula 

V~E + F = 2{l-g) (B15) 

relating the genus number g of the manifold on which the 
system is embedded to the number of local faces {F) , edges 
(E), and vertices {V) affords us with a knowledge of the 
number of independent (d = 1) symmetry operators (loops 
around independent cycles) that the system may have that 
cannot be written in terms of local operators. This and 
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